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Abstract 

Given a graph G and two vertices s and t in it, graph reachability is the problem of checking 
whether there exists a path from s to t in G. We show that reachability in directed layered 
planar graphs can be decided in polynomial time and 0(n'^) space, for any e > 0. The previous 
best known space bound for this problem with polynomial time was approximately 0{^/n) space 
|INP+13l . 

Deciding graph reachability in SC is an important open question in complexity theory and 
in this paper we make progress towards resolving this question. 


1 Introduction 

Given a graph and two vertices s and t in it, the problem of determining whether there is a path 
from s to t in the graph is known as the graph reachability problem. Graph reachability problem is 
an important question in complexity theory. Particularly in the domain of space bounded computa¬ 
tions, the reachability problem in various classes of graphs characterize the complexity of different 
complexity classes. The reachability problem in directed and undirected graphs, is complete for the 
classes non-deterministic log-space (NL) and deterministic log-space (L) respectively [LP82[|Rei08] . 
The latter follows due to a famous result by Reingold who showed that undirected reachability is in 
L [R,ein8| . Various other restrictions of reachability has been studied in the context of understand¬ 
ing the complexity of other space bounded classes (see |RTV06llGRVllllLan97| h Wigderson gave a 
fairly comprehensive survey that discusses the complexity of reachability in various computational 
models |Wig92| . 

The time complexity of directed reachability is fairly well understood. Standard graph traversal 
algorithms such as DFS and BPS solve this problem in linear time. We also have a 0(log^ n) space 
algorithm due to Savitch [SavTOj . however it requires time. The question, whether there 

exists a single algorithm that decides reachability in polynomial time and polylogarithmic space 
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is unresolved. In his survey, Wigderson asked whether it is possible to design a polynomial time 
algorithm that uses only O(n^) space, for some constant e < 1 |Wig92| . This question is also still 
open. In 1992, Barnes, Buss, Ruzzo and Schieber made some progress on this problem and gave an 
algorithm for directed reachability that requires polynomial time and 0(n/2^^°®’^) space [BBR.S92] . 

Planar graphs are a natural topological restriction of general graphs consisting of graphs that 
can be embedded on the surface of a plane such that no two edges cross. Grid graphs are a subclass 
of planar graphs, where the vertices are placed at the lattice points of a two dimensional grid and 
edges occur between a vertex and its immediate adjacent horizontal or vertical neighbor. 

Asano and Doerr provided a polynomial time algorithm to compute the shortest path (hence 
can decide reachability) in grid graphs which uses space, for any small constant e > 0 

[ADllj . Imai et al extended this to give a similar bound for reachability in planar graphs [INP'*~13] . 
Their approach was to use a space efficient method to design a separator for the planar graph 
and use divide and conquer strategy. Note that although it is known that reachability in grid 
graphs reduces to planar reachability in logspace, however since this class (polynomial time and 
Q(j.ji/2+e) gpg^gg^ ig closed Under logspace reductions, planar reachability does not follow from 
grid graph reachability. Subsequently the result of Imai et al was extended to the class of high-genus 
and H-minor-free graphs [CPT+14| . Recently Asano et al gave a 0{y/n) space and polynomial 
time algorithm for reachability in planar graphs, thus improving upon the previous space bound 
[AKNWH] . More details on known results can be found in a recent survey article [Vinl4| . 

In another line of work, Kannan et al gave a 0{n^) space and polynomial time algorithm for 
solving reachability problem in unique path graphs [KKROS] . Unique path graphs are a generaliza¬ 
tion of strongly unambiguous graphs and reachability problem in strongly unambiguous graphs is 
known to be in SC (polynomial time and polylogarithmic space) |B JLR911 ICoo79] . Reachability in 
strongly unambiguous graphs can also be decided by a 0(log^ n/log log n) space algorithm, how¬ 
ever this algorithm requires super polynomial time [AL98] . SC also contains the class randomized 
logspace or RL [Nis95] . We refer the readers to a recent survey by Allender [AllOTj to further under¬ 
stand the results on the complexity of reachability problem in UL and on certain special subclasses 
of directed graphs. 

Our Contribution 

We show that reachability in directed layered planar graphs can be decided in polynomial time and 
0{n^) space for any constant e > 0. A layered planar graph is a planar graph where the vertex set 
is partitioned into layers (say Lq to Lm) and every edge occurs between layers Lj and Lj+i only. 
Our result significantly improves upon the previous space bound due to |INP'*~13] and [AKNWl^ 
for layered planar graphs. 

Theorem 1. For every e > 0, there is a polynomial time and 0{n^) space algorithm that decides 
reachability in directed layered planar graphs. 

Reachability in layered grid graphs is in UL which is a subclass of NL |ABC^09| . Subsequently 
this result was extended to the class of all planar graphs [BTVn9] . Allender et al also gave some 
hardness results the reachability problem in certain subclasses of layered grid graphs. Specihcally 
they showed that, ILGGR is hard for NC^ and IILGGR is hard for TC'^ [ABC^Ohj . Both these 
problems are however known to be contained in L though. 

Firstly we argue that its enough to consider layered grid graphs (a subclass of general grid 
graphs). We divide a given layered grid graph into a courser grid structure along k horizontal and 
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k vertical lines (see Figure dl). We then design a modified DFS strategy that makes queries to the 
smaller graphs defined by these gridlines (we assume a solution in the smaller graphs by recursion) 
and visits every reachable vertex from a given start vertex. The modified DFS stores the highest 
visited vertex in each vertical line and the left most visited vertex in each horizontal line. We use 
this information to avoid visiting a vertex multiple number of times in our algorithm. We choose 
the number of horizontal and vertical lines to divide the graph appropriately to ensure that the 
algorithm runs in the required time and space bound. 

The rest of the paper is organized as follows. In Section [2l we give some basic definitions and 
notations that we use in this paper. We also state certain earlier results that we use in this paper. 

In Section [3l we give a proof of Theorem [TJ 

2 Preliminaries 

We will use the standard notations of graphs without defining them explicitly and follow the 
standard model of computation to discuss the complexity measures of the stated algorithms. In 
particular, we consider the computational model in which an input appears on a read-only tape 
and the output is produced on a write-only tape and we only consider an internal read-write tape 
in the measure of space complexity. Throughout this paper, by log we mean logarithm to the base 
2. We denote the set {1,2, • • • ,n} by [n]. Given a graph G, let V{G) and E{G) denote the set of 
vertices and the set of edges of G respectively. 

Definition 1 (Layered Planar Graph). A planar graph G = {V,E) is referred as layered planar if 
it is possible to represent V as a union of disjoint partitions, V = Vi U V 2 U • • • U14, for some k > 0, 
and there is a planar embedding of edges between the vertices of any two consecutive partitions Vi 
and Fj+i and there is no edge between two vertices of non-consecutive partitions. 

Now let us define the notion of layered grid graph and also note that grid graphs are by definition 
planar. 

Definition 2 (Layered Grid Graph). A directed graph G is said to be a n x n grid graph if it can 
be drawn on a square grid of size n x n and two vertices are neighbors if their Li-distance is one. 

In a grid graph a edge can have four possible directions, i.e., north, south, east and west, but if we 
are allowed to have only two directions north and east, then we call it a layered grid graph. 

We also use the following result of Allender et al to simplify our proof [ABG'*~n^ . 

Proposition 2 l [ABG'*~n9] L Reachability problem in directed layered planar graphs is log-space 
reducible to the reachability problem in layered grid graphs. 

2.1 Class nSC and its properties 

TISP(t(n), s(n)) denotes the class of languages decided by a deterministic Turing machine that runs 
in 0{t{n)) time and 0{s{n)) space. Then, SC = TISP(n‘^(^), (logn)^*^^)). Expanding the class SC, 
we define the complexity class nSC (short for near-SC) in the following definition. 

Definitions (Complexity Glass near-SC or nSC). For a fixed e > 0, we define nSC^ := TISP(n‘^(^\ n*^). 
The complexity class nSC is defined as 

nSC := Pi nSCe. 

e>0 
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We next show that nSC is closed under log-space reductions. This is an important property of 
the class nSC and will be used to prove Theorem [TJ Although the proof is quite standard, but for 
the sake of completeness we provide it here. 

Theorem 3. If A <i B and B G nSC, then A € nSC. 

Proof. Let us consider that a log-space computable function / be the reduction from A to B. It is 
clear that for any x £ A such that \x\ = n, |/(x)| < n'^, for some constant c > 0. We can think that 
after applying the reduction, f{x) appears in a separate write-once output tape and then we can 
solve /(x), which is an instance of the language B and now the input length is at most n‘^. Now 
take any e > 0 and consider e' = | > 0. B £ nSC implies that B £ nSC^/ and as a consequence, 
A £ nSCg. This completes the proof. □ 

3 Reachability in Layered Planar Graphs 

In this section we prove Theorem [TJ We show that the reachability problem in layered grid graphs, 
(denoted as LGGR) is in nSC (Theorem HJ). Then by applying Proposition [2] and Theorem |3| we 
have the proof of Theorem [TJ 

Theorem 4. LGGR £ nSC. 

To establish Theorem [TJ we define an auxiliary graph in Section 13.11 and give the required 
algorithm in Section 13.21 

3.1 The Auxiliary Graph H 

Let G be a n X n layered grid graph. We denote the vertices in G as (i, j), where 0 < i, j < n. Let 
khe a parameter that determines the number of pieces in which we divide G. We will fix the value 
of k later to optimize the time and space bounds. Assume without loss of generality that k divides 
n. Given G we construct an auxiliary graph H as described below. 

Divide G into many blocks of dimension n/k x n/k. More formally, the vertex set of H is 

V{H) := I either i or j is a non-negative multiple of n/k.} 

Note that V{H) C V{G). We consider many blocks Gi,G 2 ,--- ,Gfc 2 , where a vertex (i,j} £ 
V(Gi) if and only if < i < (i' -|- l)f and < j < (/ + l)f > for some integer i' > 0 and j' > 0 
and the vertices for which any of the four inequalities is strict will be referred as boundary vertices. 
Moreover, we have I = i' ■ k + f -|- 1. E{Gi) is the set of edges in G induced by the vertex set V(Gi). 

For every i £ [/s-l-l], let and Ly{i) denote the set of vertices, Lfi{i) = = (*~l)f} 

and Ly{i) = = {i — When it is clear from the context, we will also use L/i(i) and 

Ly{i) to refer to the corresponding gridline in H. Observe that H has k + 1 vertical gridlines and 
A; -|- 1 horizontal gridlines. 

For every pair of vertices u,v £ V{Gi) n V{H), for some I, add the edge {u,v) to E{H) if and 
only if there is a path from u to u in G/, unless u,v £ Ly{i) ov u,v £ Lh{i), for some i. Also for 
every pair of vertices u,v £ V (Gi), for some I, such that u = a-iid v = ( 12 , j 2 ), where ji = 

and j 2 = (f + l)^, for some f or ii = i'j and 12 = {f + l)^, for some i', we add edge between u 
and V in the set E[H) if and only if there is a path from rt to u in Gi and we call such vertices as 
corner vertices. 
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Figure 1: An example of layered grid graph and its decomposition into blocks 

Before proceeding further, let us introduce a few more notations that will be used later. For 
j G [k], let Lh{i,j) denote the set of vertices in Lh{i) in between L^{j) and L^{j + 1). Similarly 
we also define Ly{i,j) (see Figured]). For two vertices x,y G Ly{i), we say x y if x is below y in 
Ly{i). For two vertices x,y G we say x -< y if x is right of y in Lh{i). 


Lemma 5. Let u and v be two boundary vertices in G such that they belong to two different vertical 
or horizontal gridlines or u G Gi and v G Gj, for i j. There is a path from u to v in G if and 
only if there is path from u to v in the auxiliary graph H. 

Proof. As every edge (a, b) in H corresponds to a path from a to 6 in G, so if-part is trivial to see. 
Now for the only-if-part, consider a path P from u to x in G. P can be decomposed as P 1 P 2 ■ ■ ■ Pr, 
such that Pi is a path from x, G V{Gi) to Xj+i G V{Gi), where it re-enters V{Gi) for the next time, 
for some I, and is of following two types: 

1. Xi and Xj+i belong to different horizontal or vertical gridlines; or 

2. Xi and Xj+i are two corner vertices. 

Now by the construction H, for every i, there must be an edge (xj,Xi+i) in H for both the 
above cases and hence there is a path from x to x in Ff as well. □ 

Now we consider the case when u and x belong to the same vertical or horizontal gridlines. 

Claim 6. Let u and v be two vertices contained in either Ly{i) or Lh(i) for some i. Then deciding 
reachability between u and v in G can be done in logspace. 
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Proof. Let us consider that u,v G Ly{i), for some i. As the graph G under consideration is a layered 
grid graph, if there is a path between u and v, then it must pass through all the vertices in L^{i) 
that lies in between u and v. Hence just by exploring the path starting from u through Ly(i), we 
can check the reachability and it is easy to see that this can be done in logspace, because the only 
thing we need to remember is the current vertex in the path. Same argument will also work when 
u,v G Lh{i), for some i and this completes the proof. □ 

Now we argue on the upper bound of the length of any path in the auxiliary graph H. 

Lemma 7. For any two vertices u,v G V{H), any path between them is of length at most 2k + 1. 

Proof. Consider any two vertices u,v G V (H) and a path u = xiX 2 ■ ■ ■ Xr = v, from u to v, denoted 
as P. Now let us consider a bipartite (undirected) graph K, where V{K) = Au B such that 
A = {xi\i G [r]} and B = {L^(i), L/i(j)|i, j G [k + 1]}. We add an edge (a, 6) in E{K) if and 
only if a G 6, where b = L.„{i) or b = Lh{i), for some i. Now observe that by the constrnction 
of H, each Xj belongs to different Ly{i) or Lfi{i) unless xj is some corner vertex and in that case 
Xj G Ly{i),Lh{i'), for some i and i'. Moreover, if Xj G Ly{i) (or Ly{i)), but is not a corner vertex, 
then Xj+i cannot be a corner vertex on Ly(i) (or Lh{i)). As a conseqnence for every snbset Sa F A, 
its neighbor set N{Sa) ■= {b G B\3a G Sa snch that (a, 6) G E{K)} C B satisfies the condition 
that |A'(S'^)| > |5^|. Now we apply the Hall’s Theorem |LP86j . which states that a bipartite 
(undirected) graph G = {AU B,E) has a matching if and only if for every S' C A, |A^(S)| > |S|. 
Hence there is a matching between A and B and as \B\ < 2{k + 1), so is |A|. This shows that the 
path P is of length at most 2k + 1. □ 

3.2 Description of the Algorithm 

We next give a modified version of DFS that starting at a given vertex, visits the set of vertices 
reachable from that vertex in the graph H. At every vertex, the traversal visits the set of outgoing 
edges from that vertex in an anticlockwise order. 

In our algorithm we maintain two arrays of size A: + 1 each, say Ay and A/^, one for vertical and 
the other for horizontal gridlines respectively. For every i G [A: + 1], Ay{i) is the topmost visited 
vertex in Ly{i) and analogously Ah{i) is the leftmost visited vertex in Lh{i). This choice is guided 
by the choice of traversal of our algorithm. More precisely, we cycle through the outgoing edges of 
a vertex in an anticlockwise order. 

We perform a standard DFS-like procedure, using the tape space to simulate a stack, say S. 
S keeps track of the path taken to the current vertex from the starting vertex. By Lemma [7l the 
maximum length of a path in H is at most 2A: + 1. Whenever we visit a vertex in a vertical gridline 
(say Ly{i)), we check whether the vertex is lower than the f-th entry of Ay. If so, we return to 
the parent vertex and continue with its next child. Otherwise, we update the f-th entry of Ay to 
be the current vertex and proceed forward. Similarly when visit a horizontal gridline (say Lh{i)), 
we check whether the current vertex is to the left of the f-th entry of A^. If so, we return to the 
parent vertex and continue with its next child. Otherwise, we update the f-th entry of A^ to be 
the current vertex and proceed. The reason for doing this is to avoid revisiting the subtree rooted 
at the node of an already visited vertex. The algorithm is formally defined in Algorithm [TJ 

Lemma 8. Let Gi be some block and let x and y be two vertices on the boundary of Gi such that 
there is a path from x to y in G. Let x' and y' be two other boundary vertices in Gi such that (i) 


6 




Input : The auxiliary graph H, two vertices u,v G V{H) 

Output : YES if there is a path from u to v] otherwise NO 

1 Initialize two arrays Ay and Ah and a stack 5; 

2 Initialize three variables curr, prev and next to NULL; 

3 Push u onto S] 

4 while S is not empty do 

5 curr <r- top element of S'; 

6 next G- neighbor of curr next to prev in counter-clockwise order; 

7 while next / NULL do 

/* cycles through neighbors of curr */ 

8 if next = V then 

9 I return YES; 

10 end 

11 if next G Ly{i) for some i and Ay\i\ -< next then 

12 Ay[i] G- next] 

13 break; 

14 end 

15 if next G Lh{i) for some i and Ah{i\ -< next then 

16 Ah[i] G- next] 

17 break; 

18 end 

19 prev G- next] 

20 next G- neighbor of curr next to prev in counter-clockwise order; 

/* NULL if no more neighbors are present */ 

21 end 

22 if next = NULL then 

23 remove curr from S'; 

24 prev G- curr] 

25 else 

26 add next to S] 

27 prev G- NULL; 

28 end 

29 end 

30 return NO; 

Algorithm 1: AlgoLGGR: Algorithm for Reachability in the Auxiliary Graph H 
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there is a path from x' to y' in G and (ii) x' lies on one segment of the boundary of Gi between 
vertiees x and y and y' lies on the other segment of the boundary. Then there is a path in G from 
X to y' and from x' to y. Hence, if {x,y) and {x',y') are present in E{H) then so are {x,y') and 
{x',y)- 

Proof. Since G is a layered grid graph hence the paths x to y and x' to y' must lie inside Gi. 
Also because of planarity, the paths must intersect at some vertex in Gi. Now using this point of 
intersection, we can easily show the existence of paths from x to y' and from x' to y. □ 

Lemma [9] will prove the correctness of Algorithm [TJ 

Lemma 9. Let u and v be two vertices in H. Then starting at u Algorithmic visits v if and only 
if V is reachable from u. 

Proof. It is easy to see that every vertex visited by the algorithm is reachable from u since the 
algorithm proceeds along the edges of H. 

By induction on the shortest path length to a vertex, we will show that if a vertex is reachable 
from u then the algorithm visits that vertex. Let B^{u) be the set of vertices reachable from u 
that are at a distance d from u. Assume that the algorithm visits every vertex in Let x 

be a vertex in Bd{u). Without loss of generality assume that x is in L.„{i,j) for some i and j. A 
similar argument can be given if x belongs to a horizontal gridline. Further, let x lie on the right 
boundary of a block Gi. Let Wx = {tc E B(i-i{u)\{w,x) E E{H)}. Note that by the definition of 
H, all vertices in Wx lie on the bottom boundary or on the left boundary of Gi. 

Suppose the algorithm does not visit x. Since x is reachable from u via a path of length d, 
therefore Wx is non empty. Let w be the first vertex added to Wx by the algorithm. Then w is 
either in L^ij), or in — 1). Without loss of generality assume w is in Lh{j). Let 2 be the 
value in Ay{i) at this stage of the algorithm (that is when w is the current vertex). Since x is not 
visited hence x < z. Also this implies that z was visited by the algorithm at an earlier stage of the 
algorithm. Let w' be the ancestor of z in the DFS tree such that w' is in Lh{j). There must exist 
such a vertex because z is above the j-th horizontal gridline, that is Lh{j). 

Suppose if w' lies to the left of w then by the description of the algorithm, w is visited before 
w'. Hence x is visited before z. On the other hand, suppose if w' lies to the right of w. Clearly 
w' cannot lie to the right of vertical gridline Lx{i) since z is reachable from w' and z is in Lx{i). 
Let w" be the vertex in Lh{j + 1) such that w" lies in the tree path between w' and z (See Figure 
[2]). Observe that all four vertices lie on the boundary of Gi. Now by applying Lemma [8] to the 
four vertices w, x, w' and w” we conclude that there exists a path from w' to x as well. Since 
x < z, X must have been visited before z from the vertex w'. In both cases, we see that z cannot 
be Ay(i) when w is the current vertex. Since z was an arbitrary vertex such that x < z, the lemma 
follows. □ 


We next show Lemma fTOl which will help us to achieve a polynomial bound on the running time 
of Algorithm [TJ 

Lemma 10. Every vertex in the graph H is added to the set S at most once in Algorithmic 


Lh{j + 1 ) 



LhU) 

L,{i - 1) ^ L,{i) 

Figure 2: Crossing between two paths 


Proof. Observe that a vertex u in Ly{i) is added to S only if Ay{i) -< u, and once u is added, Ay{i) 
is set to u. Also during subsequent stages of the algorithm, if Ay{i) is set to v, then u P v. Hence 
u -< Ay{i). Therefore, u cannot be added to S again. 

We give a similar argument if u is in Lh{i). Suppose if u is in Ly{i) for some i and Lh{j) for 
some j, then we add u only once to S. This check is done in Line 16 of Algorithm [TJ However we 
update both Ay{i) and Ah{j). □ 

Algorithm [1] does not explicitly compute and store the graph H. Whenever it is queried with 
a pair of vertices to check if it forms an edge, it recursively runs a reachability query in the 
corresponding block and produces an answer. The base case is when a query is made to a grid of 
size k X k. For the base case, we run a standard DFS procedure on the k x k size graph. 

In every iteration of the outer while loop (Lines 4 - 29) of Algorithm [H either an element is 
added or an element is removed from S. Hence by Lemma flOl the loop iterates at most Ank times. 
The inner while loop (Lines 7 - 21), cycles through all the neighbors of a vertex and hence iterates 
for at most 2n/k times. Each iteration of the inner while loop makes a constant number of calls to 
check the presence of an edge m. a n/k x n/k sized grid. Let 'T{n) and S{n) be the time and space 
required to decide reachability in a layered grid graph of size n x n respectively. Then, 


T{n) = 


8n'^{T{n/k) + 0(1)) if n > k 
0{k‘^) otherwise. 


Hence, 


T{n) = O (n '°e 


log n 

k 


Since we do not store any query made to the smaller grids, therefore the space required to check 
the presence of an edge in H can be reused. Ay and A^ are arrays of size k + 1 each. By Lemma 
0 the number of elements in S at any stage of the algorithm is bounded by 2/c + 1. Therefore, 

\S{n/k) + 0{k\ogn) li n > k 
S{n) = < 

I 0{k‘^) otherwise. 


Hence, 


5(n) = O 


log k 


log^ n + k^ 
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Now given any constant e > 0, if we set k = then we get T{n) = and S{n) = 0{n^). 

This proves Theorem [H 
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